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Abstract 

We have investigated the electrical transport properties of Dirac electrons in a monolayer 
graphene sheet in the presence of a perpendicular magnetic field that is modulated weakly and 
periodically along one direction. We find that the Landau levels broaden into bands and their 
width oscillates as a function of the band index and the magnetic fieldWe determine the a yy 
component of the magnetoconductivity tensor for this system which is shown to exhibit Weiss 
oscillations. We also determine analytically the asymptotic expressions for a yy .We compare these 
results with recently obtained results for electrically modulated graphene as well as those for mag- 
netically modulated conventional two-dimensional electron gas (2DEG) system. We find that in the 
magnetically modulated graphene system cosidered in this work, Weiss oscillations in a yy have a 
reduced amplitude compared to the 2DEG but are less damped by temperature while they have 
a higher amplitude than in the electrically modulated graphene system. We also find that these 
oscillations are out of phase by tt with those of the electrically modulated system while they are 
in phase with those in the 2DEG system. 



1 



I. INTRODUCTION 



The successful preparation of monolayer graphene has allowed the possibility of study- 
ing the properties of electrons in graphene The nature of quasiparticles called Dirac 
electrons in these two-dimensional systems is very different from those of the conventional 
two-dimensional electron gas (2DEG) realized in semiconductor heterostructures. Graphene 
has a honeycomb lattice of carbon atoms. The quasiparticles in graphene have a band 
structure in which electron and hole bands touch at two points in the Brillouin zone. At 
these Dirac points the quasiparticles obey the massless Dirac equation. In other words, they 
behave as massless Dirac particles leading to a linear dispersion relation = vk ( with the 
characteristic velocity v ~ 10 6 m/s). This difference in the nature of the quasiparticles in 
graphene from conventional 2DEG has given rise to a host of new and unusual phenomena 
such as anamolous quantum Hall effects and a it Berry phase Earlier it was found 
that if conventional 2DEG is subjected to artificially created periodic potentials in the sub- 
micrometer range it leads to the appearence of Weiss oscillations in the magnetoresistance. 
This type of electrical modulation of the 2D system can be carried out by depositing an array 



of parallel metallic strips on the surface or through two interfering laser beams 

Besides the fundamental interest in understanding the electronic properties of graphene 
there is also serious suggestions that it can serve as the building block for nanoelectronic 
devices [6fl. Since Dirac electrons can not be confined by electrostatic potentials due to the 

n 

Klein's paradox it was suggested that magnetic confinement be considered |7j]. Technology 
for this already exists as the required magnetic field can be created by having ferromagnetic 
or superconducting layers beneath the substrate s|. 

In conventional 2DEG systems, electron transport in the presence of magnetic barri- 
ers and superlattices has continued to be an active area of research js]. Recently, electrical 
transport in graphene in the presence of electrical modulation was considered and theoretical 
predictions made 111] .Along the same lines, in this work we investigate low temperature mag- 
netotransport of Dirac electrons in a single graphene layer subjected to a one-dimensional 
(ID) magnetic modulation. The perpendicular magnetic field is modulated weakly and 
periodically along one direction 
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II. FORMULATION AND ENERGY SPECTRUM 



We consider two-dimensional Dirac electrons in graphene moving in the x-y-plane. The 
magnetic filed (B) is applied along the z-direction perpendicular to the graphene plane. The 
perpendicular magnetic field B is modulated weakly and periodically along one direction 
such that B = {B + Bqcos(Kx))^. Here Bq is the strength of the magnetic modulation. 
In this work we consider the modulation to be weak such that B << B. We consider the 
graphene layer within the single electron approximation. The low energy excitations are 



described by the two-dimensional (2D) Dirac like Hamiltonian (h = c = 1 here) 



QflQ 



H = v*o>.{-iV +eA). (1) 

Here ^a* = { a x , are the Pauli matrices and v characterizes the electron velocity. We 
employ the Landau gauge and write the vector potential as A = (0, Bx + (B /K) sin^x), 0) 
where K = 27r/a and a is the period of the modulation. The Hamiltonian given by Eq. (1) 
can be expressed as 

y 

H = — iv a . V + ev It^Bx + ev tx^ — sin(i^:r). (2) 

K 

The above Hamiltonian can be written as 

H = H + H[ (3) 
where H is the unmodulated Hamiltonian given as 

Hq = —iv^.'V + evVyBx 

and 



H' = eva v — sm(Kx). 
K 

The Landau level energy eigenvalues without modulation are given by 

e(n) = uj g ^/n (4) 



where n is an integer and u g = v\/2eB. As has been pointed out Hj the Landau level 
spectrum for Dirac electrons is significantly different from the spectrum for electrons in 
conventional 2DEG which is given as e(n) = u c (n + 1/2), where u c = eB/m is the cyclotron 
frequency. 



The eigenfunctions without modulation are given by 



e 



ik y y 



( 



where 



-i^n^Kx + X )/l] 



y $ n [(x + x )/i] 

x 2 /2 



(5) 



= 6 H n (x) (6) 



where / = ^Jl/eB is the magnetic length, xq = l 2 k y , L y is the y-dimension of the graphene 
layer and H n (x) are the Hermite polynomials. Since we are considering weak modulation 
Bq << B, we can apply standard perturbation theory to determine the first order corrections 
to the unmodulated energy eigenvalues in the presence of modulation 

OO Ly 

Ae^ = JdxJ dyK, ky (r)H>(x)y n>ky (r) (7) 

-oo 

with the result 

Ae^ ky = u cos(Kx ) (2 v ^ e -"/ 2 [L n _ 1 ( M ) - L n (u)}) (8) 

where co = j^zf, u = K 2 l 2 /2 and L n {u) are the Laguerre polynomials. Hence the energy 
eigenvalues in the presence of modulation are 

e(n, k y ) = e(n) + Ae n ky = co g y/n + u cos(Kx )G n (9) 

with G n (u) = 2 y /ne~ u / 2 [L n _ 1 (u) — L n (u)]. We observe that the degeneracy of the Landau 
level spectrum of the unmodulated system with respect to k y is lifted in the presence of 
modulation with the explicit presence of k y in xq. The n = landau level is different from 
the rest as the energy of this level is zero and electrons in this level do not contribute to 
diffusive conductivity calculated in the next section. The rest of the Landau levels broaden 
into bands. The Landau bandwidths ~ G n oscillates as a function of n since L n (u) are 
oscillatory functions of the index n. 

Before we begin the calculation of electrical conductivity it is necessary to discuss the 
regime of validity of the perturbation theory presented above. For large n the level spacing 
given by Eq. (4) goes as uj g (^fn — y/(n — 1)) — > ^g^/K anc ^ ^ ne width of the nth level 
given by Eq.(8) goes as 2cu n 1 ^ 2 . There is therefore a value of n at which the width becomes 
equal to the spacing and the perturbation theory is no longer valid. This occurs when 
n max = v^tt 2 -^ where B' — — 0.0054T for a = 350nm. For a fixed electron density and 



the period of modulation this suggests the maximum value of the magnetic modulation Bo 
above which it is necessary to carry out a more sophisticated analysis. 



III. ELECTRICAL CONDUCTIVITY WITH PERIODIC MAGNETIC MODULA- 
TION 

To calculate the electrical conductivity in the presence of weak magnetic modulation we 
use Kubo formula to calculate the linear response to an applied external field. In a magnetic 
field, the main contribution to Weiss oscillations comes from the scattering induced migration 
of the Larmor circle center. This is diffusive conductivity and we shall determine it following 
the approach in 10|, QJ| where it was shown that the diagonal component of conductivity 
a yy can be calculated by the following expression in the case of quasielastic scattering of 
electrons 



°yy = T-rY.f^ 1 ~ f( E <MEM) 2 (io) 



L x , L y , are the dimensions of the layer, (3 = is the inverse temperature with ks the 
Boltzmann constant, f(E) is the Fermi Dirac distribution function, t(E) is the electron 
relaxation time and £ denotes the quantum numbers of the electron eigenstate. The diagonal 
component of the conductivity a yy is due to modulation induced broadening of Landau bands 
and hence it carries the effects of modulation in which we are primarily interested in this 
work. a xx does not contribute as the component of velocity in the x-direction is zero here. 
The collisional contribution due to impurities is not taken into account in this work. 
The summation in Eq.(10) over the quantum numbers ( can be written as 



x V C g n=0 n=0 

The component of velocity required in Eq.(10) can be calculated from the following expres- 
sion 

d 



v y = W y e{n ' ky) - (12) 

Substituting the expression for e(n, k y ) obtained in Eq.(9) into Eq.(12) yields 

v y = sm(Kxo)G n (u) (13) 



With the results obtained in Eqs.(ll), (12) and (13) we can express the diffusive contribution 
to the conductivity given by Eq.(lO) as 

<Tyy = A <l> (14) 

where 

A = 2c^e 2 r/3 (15) 
and the dimensionless conductivity $ is given as 

$ = Aue ~ u it i^fi)p [Ln - x{u) - Ln{u)]2 (i6) 

where g(E) = exp[/3(E — Ep\ and Ep is the Fermi energy. 

In Fig.(l), we plot the dimensionless conductivity given by Eq.(16) as a function of 
inverse magnetic field at temperature T = 6K and electron density n e = 3 x 10 u cm~ 2 . The 
dimensionless magnetic field is introduced given as b = % with B' = -V In the region of 
high magnetic field we can see SdH oscillations superimposed on the Weiss oscillations. 



IV. ASYMPTOTIC EXPRESSIONS 



To get a better understanding of the results of the previous section we will consider the 
asymptotic expression of conductivity where analytic results in terms of elementary functions 



can be obtained following ll|]. We shall compare the asymptotic results for the dimensionless 
conductivity obtained in this section with the results obtained for a magnetically modulated 
conventional 2DEG system. We shall also compare these results with those of graphene that 
is subjected to only the electric modulation 

The asymptotic expression of dimensionless conductivity can be obtained by using the 
following asymptotic expression for the Laguerre polynomials 

exp-" /2 L n (u) -> 1 cos(2v/mZ - E). (17) 

Note that the asymptotic results are valid when many Landau Levels are filled. We now 
take the continuum limit: 



i (ieV ^ /r 2 



n -->- 2 {-) ■s— > W l EdE (18) 



o 



to express the dimensionless conductivity in Eq.(16) as the following integral 



u ( I 



oo 

/ dE [g{E) + l)Y sin2 (V2V^) sm 2 (2 



7T . 



nu — 



(19) 



where u = 2n 2 /b. 

Now assuming that the temperature is low such that (5 <C Ep and replacing E 
E F + we rewrite the above integral as 



5P a . 2 / 7T 

sin ' 
1 p 



<ise 



(e s + 1) 



2 27rp 7r 27ra 
sin (— 1 — —s) 



4 vbp 



(20) 



where p = = kpa = ^/2ixn e a is the dimensionless Fermi momentum of the electron. To 
obtain an analytic solution we have also replaced E by Ep in the above integral except in 
the sine term in the integrand. 

The above expression can be expressed as 



W a . 2 / 7T 
sin — 



vb 2 (3 



ds 



pj J coslr(s/2) 



2 27rp 7r 27ra 
sin — H — —s] 
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The above integration can be performed by using the following identity 



12|: 



(21) 



dx 



cos ax 



an 



cosh 2 /^ 2/5 2 sinh(a7r/2/3) 



with the result 
$ = 

where ksTp, = 



2p 2 T . 2 /vr 

2^ S111 - 



1 - A, £) +2A (^ lsini 
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r = 4zrfa and = _x _(*— >oo) _ >= 2a;e -s. 



(22) 



(23) 



V. COMPARISON WITH MAGNETICALLY MODULATED CONVENTIONAL 
TWO-DIMENSIONAL ELECTRON GAS SYSTEM 



We start by comparing the energy spectrum and velocity expression obtained in Eq.(9) 
and Eq.(13) with similar expressions for the conventional 2DEG where the electron spec- 
trum is parabolic [9j. For the energy spectrum, we find that the Landau level spectrum is 
significantly different from that of standard electrons in conventional 2DEG. The first term 



LUgy/n with ujg = v\ / 2eB in Eq.(9) has to be compared with u c (n + 1/2) with uj c = eB/m 
for standard electrons. The modulation effects are in the second term where the essential 
difference is in the structure of the function G n {u) = 2yfne~ u l 2 [L n -i{u) — L n (u)). We find 
that there are essentially two basic differences: Firstly, for Dirac electrons we have a dif- 
ference of two successive Laguerre polynomials whereas we had the sum of the Laguerre 
polynomials in the corresponding term for standard electrons in 2DEG. Secondly, the ex- 
pression for Dirac electrons is multiplied by the square root of the Landau band index ^Jn 
that was absent in the expression for standard electrons. The above mentioned differences 
in the G n {u) function cause the velocity expression for the Dirac electrons given by Eq.(13) 
to be different from that of the standard electrons. 

As expected, these differences in the energy spectra and velocities lead to different results 
for the diffusive conductivity in the two cases. We now compare the results obtained for the 
asymptotic expression of the diffusive conductivity a yy . To make this comparison possible 
we first express $ given by Eq.(23) as 

* = dhlTsin 2 - )F (24) 



ve 2 B 2 a \p 

where F — 1 — A (^gj + sm2 [2vr (| — |)] . We now compare the results for di- 

mensionless conductivity obtained in Eq.(24) with those presented in Eq. (18) of {^(a). We 
find that the result for graphene (Dirac electrons) system differs from that of the conven- 
tional 2DEG (standard electrons) system by a factor ^ sin 2 . For the system under 

consideration p ~ 50 and in this limit if we take sin 2 y~j — > (^j it yields the factor |. 
Hence we conclude that the amplitude of the oscillations in the conductivity will be reduced 
by this factor in the magnetically modulated graphene system compared to the conventional 
2DEG under the same conditions. For the parameters considered in this work, the conduc- 
tivity is larger by a factor of ~ 44 (at magnetic field 0.5T) in the 2DEG system compared 
to graphene. In Fig. (2) we plot the dimensionless conductivity versus inverse magnetic field 
for magnetically modulated graphene and 2DEG. Note that in Fig. (2) the dimensionless 
conductivity for conventional 2DEG is rescaled by a factor .023. 

The temperature scale for damping of Weiss oscillations in graphene can be obtained from 
Eq.(23) and is characterized by To given above while the characteristic tempererature for 
2DEG is given in [9j(a) as k B T a = (^cj c /47t 2 ) ak F . Comparing T D and T a we obtain ^ = ^ 
where Vp is the Fermi velocity in 2DEG. This shows that Weiss oscillations in graphene are 



less damped with temperature compared to 2DEG due to the difference in Fermi velocities 
in the two systems. 



VI. COMPARISON WITH ELECTRICALLY MODULATED GRAPHENE SYS- 
TEM 



We will now compare the results obtained in this work with results obtained in [11] 
for the case of electrically modulated graphene system. We will first compare the energy 
spectrum in the two cases. The difference in the energy spectrum due to modulation effects 
was obtained in Eq.(8). If we compare this result with the corresponding expression for 
the electrically modulated case, we find the following differences: Firstly, in the magnetic 
modulation case we have a difference of two successive Laguerre polynomials whereas we 
had the average of two successive Laguerre polynomials in the electric case. Secondly, in 
the magnetic modulation case the energy eigenvalues are multiplied by the square root of 
the Landau band index y/n that was absent in the expression for the electric case. These 
differences cause the velocity expression for the Dirac electrons given by Eq.(13) to be 
different from that of electrons in the electrically modulated system. 

We now compare the expressions for dimensionless conductivity $ given by Eq. (23) with 
the electrically modulated case (Eq.(22) in [11]). We find that in the magnetically modulated 
case we have sin 2 x functions in place of cos 2 x functions for the electric case which results 
in oscillations being out of phase in the two cases. We also find that the amplitude of the 
oscillations in the magnetic case are larger by a factor of compared to the electrically 
modulated case. For the parameters considered in this work, the conductivity is larger by 
a factor of ~ 1.2 (at magnetic field 0.5T) in the magnetically modulated graphene system 
compared to the electrically modulated one. 

Exact expression of the dimensionless conductivity $ for the electric and magnetic mod- 
ulated graphene system is shown in Fig. (3) as a function of the inverse magnetic field at 
temperature T = 6K , electron density n e = 3 x 10 n cm~ 2 and period of modulation 
a = 350nm. In Fig. (3) we can clearly see that the Weiss oscillations in the dimensionless 
conductivity are enhanced, they have a larger amplitude, in the magnetically modulated 
case compared to the electrically modulated case for the same parameter values. Further- 
more, we note that the oscillations in the magnetic and electric modulated cases have a tt 
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phase shift. We also observe that in the region of high magnetic field SdH oscillations are 
superimposed on the Weiss oscillations. The oscillations are periodic in 1/5 and the period 
depends on electron density as ^JrT e in both the magnetic and electric modulated cases. The 
characteristic damping temperature is the same for both the systems. 

To better understand the increase in amplitude of Weiss oscillations in the magnetically 
modulated graphene system compared to the electrically modulated one we consider the 
difference in bandwidths in the two cases [kJ. Important feature is the additional yfn factor 
in the perturbed energy eigenvalues for the magnetically modulated case which is absent 
in the electrically modulated case. The result is that the bandwidth in the magnetically 
modulated case is approximately greater by a factor of Ay/u compared to the electrically 
modulated graphene system. 

VII. CONCLUSIONS 

In this work we have investigated the electrical transport properties of Dirac electrons 
in a monolayer graphene sheet in the presence of a perpendicular magnetic field that is 
modulated weakly and periodically along one direction. Our primary focus has been the 
study of Weiss oscillations in the diffusive magnetoconductivity a yy of this system. We have 
compared the results obtained with those obtained for magnetically modulated conventional 
2DEG system and with those of the graphene system subjected to only the electric modula- 
tion. We find that in the magnetically modulated graphene system Weiss oscillations in the 
magnetoconductivity have a reduced amplitude compared to the conventional 2DEG but are 
more robust with respect to temperature. In comparison with the electrically modulated 
graphene case, we find that the conductivity is larger in amplitude. We also find that the 
oscillations in the magnetoconductivity in graphene are 7r phase shifted with respect to the 
electrically modulated case whereas they are in phase with the conventional 2DEG subjected 
to magnetic modulation. 

VIII. ACKNOWLEDGEMENTS 

One of us (K.S.) would like to acknowledge the support of the Pakistan Science Foun- 
dation (PSF) through project No. C-QU/Phys (129). M. T. would like to acknowledge the 



10 



support of the Pakistan Higher Education Commission (HEC). 

♦Present address: Department of Physics, Blackett Laboratory, Imperial College London, 
London SW7 2AZ, United Kingdom. 



[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. 

V. Dubonos, and A. A. Firsov, Nature 438, 197 (2005); Y. Zhang, Y. W. Tan, H. L. Stormer, 

and P. Kim, 438, 201 (2005). 
[2] Y. Zheng and T. Ando, Phys. Rev. B 65, 245420 (2002); V. P. Gusynin and S. G. Sharapov, 

Phys. Rev. Lett. 95, 146801 (2005); N. M. R. Perez F. Guinea, and A. H. Castro Neto, Phys. 

Rev. B 73, 125411 (2006); M. I. Katsnelson , K. S. Novoselov, and A. K. Geim, Nat. Phys. 

2, 620 (2006); K. S. Novoselov, E. McCann, S. V. Morozov, V. I. Fal'ko, M. I. Katsnelson, U. 

Zeitler, D. Jiang, F. Schedin, and A. K. Geim, Nat. Phys. 2, 177 (2006). 
[3] D. Weiss, K. v. Klitzing, K. Ploog, and G. Weimann, Europhys. Lett., 8, 179 (1989) 
[4] R. W. Winkler, J. P. Kotthaus, and K. Ploog, Phys. Rev. Lett. 62, 1177 (1989). 
[5] R. R. Gerhardts, D. Weiss, and K. v. Klitzing, Phys. Rev. Lett. 62, 1173 (1989). 
[6] C. Berger,et.a/, Science 312, 1191 (2006). 

[7] A.De Martino, L. Dell' Anna, R. Egger, Phys. Rev. Lett. 98, 066802 (2007). 
[8] S. J. Lee et.al, Phys.Rep. 394, 1 (2004). 

[9] (a) F. M. Peeters and P. Vasilopoulos, Phys. Rev. B 47, 1466 (1993); (b) J. Shi and F. 
M. Peeters, K. W. Edmonds and B. L. Gallagher, Phys. Rev. B 66, 035328 (2002); (c) P. 
Vasilopoulos, F. M. Peeters, Superlattices and Microstructures 7, 393 (1990); (d) F. M. Peeters 
and A. Matulis, Phys. Rev. B 48, 15166 (1993); (e) D. P. Xue and G. Xiao, Phys. Rev. B 45, 
5986 (1992); (f) P. De Ye, D. Weiss, R. R. Gerhardts, M. Seeger, K. v. Klitzing, K. Eberl and 
H. Nickel, Phys. Rev. Lett. 74, 3013 (1995); (g) J. H. Ho, Y. H. Lai, Y. H. Chui and M. F. 
Lin, Nanotechnology 19, 035712 (2008). 

[10] F. M. Peeters and P. Vasilopoulos, Phys. Rev. B 46, 4667 (1992); M. Charbonneau, K. M. 
Van Vliet and P. Vasilopoulos, J. Math. Phys. 23, 318 (1982). 

[11] A. Matulis and F. M. Peeters, Phys. Rev. B 75, 125429 (2007). 

[12] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products (Academic Press, 
New York, 1980). 

11 



[13] M. Tahir, K.Sabeeh and A. MacKinnon, J. Phys.:Condens. Matter 19, 406226 (2007). 



12 




0.00 0.01 0.02 0.03 0.04 0.05 

B'/B 

Fig .(1 ) : The dimensionless conductivity versus inverse magnetic field. 



